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NMR and magnetization measurements in Li2VOSiC>4 and I^VOGeCU are reported. The analysis 
of the susceptibility shows that both compounds are two-dimensional S = 1/2 Heisenberg antifer- 
romagnets on a square lattice with a sizeable frustration induced by the competition between the 
superexchange couplings Ji along the sides of the square and J2 along the diagonal. Li2VOSi04 
undergoes a low-temperature phase transition to a collinear order, as theoretically predicted for 
J2/J1 > 0.5. Just above the magnetic transition the degeneracy between the two collinear ground 
states is lifted by the onset of a structural distortion. 



PACS numbers: 76.60.Es, 75.40. Gb, 75.10.Jm 

In recent years one has witnessed an extensive inves- 
tigation of quantum phase transition in low-dimensional 
5=1/2 Heisenberg antifcrromagnets (QHAF) as a func- 
tion of doping, magnetic field and disorder |fj . For ex- 
ample, two-dimensional QHAF (2DQHAF) have been 
widely studied to show the occurrence of a phase transi- 
tion from the renormalized classical to the quantum dis- 
ordered regime upon charge doping S . Another possibil- 
ity to drive quantum phase transitions in a 2DQHAF is to 
induce a sizeable frustration. In particular, for a square 
lattice with an exchange coupling along the diagonal J2 
about half of the one along the sides of the square J\ (see 
Fig. la), a crossover to a quantum disordered phase with 
a finite gap between the singlet ground state and the first 
excited state is expected |3]-|j|. For J2/J1 -C 0.5 a Neel 
order is envisaged, while for J%j J\ 3> 0.5 a collinear or- 
der should develop. The collinear order (see Fig. la), 
which can be considered as formed by two interpene- 
trating Neel sublattices with staggered magnetization ni 
and ri2, is characterized by an Ising order parameter 
a = ni.ri2 = ±1 ||. The two values of a correspond 
to the two collinear configurations, one with spins ferro- 
magnetically aligned along the x axis, with a magnetic 
wave- vector Q = (0, n/a), the other with spins ferromag- 
netically aligned along the y axis (Q = (7r/a,0)). At a 
certain temperature an Ising phase transition occurs and 
the system choses among the x or y collinear configura- 
tions. The precise boundaries of the Ji / J\ phase dia- 
gram for a frustrated 2DQHAF are unknown and could 
be modified by the presence of a finite third neighbour 
coupling H . These theoretical predictions have not found 
an experimental support so far, mainly due to the ab- 
sence of a system which can be regarded as a prototype 
of a frustrated 2D QHAF. 

In this letter we present NMR and magnetization mea- 
surements that prove that the isostructural compounds 
Li 2 VOSi0 4 (LSVO for short) and Li 2 VOGe0 4 (LGVO) 



||, formed by layers of V 4+ (S — 1/2) ions on a 
square lattice (see Fig. lb), are prototypes of frustrated 
2DQHAF with significant coupling between both first 
(Ji) and second (Ji) neighours. Moreover we show that 
LSVO undergoes a phase transition to a low temperature 
collinear order, as expected for J2/J1 > 0.5. The phase 
transition is triggered by a lattice distortion which lifts 
the degeneracy between the two possible collinear ground 
states and could belong to the Ising universality class. 

29 Si NMR and magnetization measurements have been 
performed on powder samples while 7 Li NMR measure- 
ments, thanks to 7 Li sensitivity, have been carried out 
also on a ~ 1 x 1 x 0.2 mm 3 LSVO single crystal. NMR 
spectra and nuclear spin-lattice relaxation rate 1/T± have 
been measured by using standard pulse sequences. The 
field-cooled magnetization M was measured with a com- 
mercial Quantum Design MPMS-XL7 SQUID magne- 
tometer. 

The structure of V A+ layers || suggests that both 
the couplings between first and second neighbours can 
be significant. It is however difficult a priori to decide 
which one should dominate: first neighbours are con- 
nected by two superexchange channels, but they are lo- 
cated in pyramids looking in opposite directions and are 
not exactly in the same plane, whereas second neighbours 
are connected by only one channel, but they are located 
in pyramids looking in the same directions and are in the 
same plane. It would thus be highly desirable to extract 
information on the relative value of these exchange in- 
tegrals from the susceptibility (\ — M/H) (see Fig. 2). 
Although the temperature dependence of the suscepti- 
bility of the J\ — Ji model is not known accurately as 
a function of J\ and Ji, it turned out to be possible to 
obtain useful information from the following considera- 
tions. If the system was not frustrated, i.e. if J\ ^> J2 
or J2 3> J\ , the susceptibility would be that of a regular 
Heisenberg AF on the square lattice with coupling J. In 
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that case, Quantum Monte Carlo simulations can be used 
to determine the temperature dependence of the suscep- 
tibility 0, and the maximum occurs at T max ~ 0.935 J. 
Since in that case the Curie- Weiss temperature = J 
one has a ratio Tmax/© — 0.935. Now, while T max is 
known very accurately from our measurements, the pre- 
cise determination of O is more problematic. A simple 
fit with X {T) = Xvv + C/{T + 0), with xvv Van-Vleck 
susceptibility, is not good enough because there is an 
important dependence of the results on the lowest tem- 
perature used in the fit. To overcome this problem, we 
have performed a fit of the high temperature part of the 
susceptibility up to third order. The coefficients are then 
consistent with a J\ — J2 model only in a small win- 
dow for the lowest temperature. Within this window, 
O depends very weakly on the lowest temperature and 
a precise estimate of the Curie- Weiss temperature can 
be achieved. The results are ~ 7AK for LSVO and 
~ 5.2K for LGVO §. Accordingly, the ratios T max /0 
are equal to 0.72 and 0.67, respectively. In both cases, 
this ratio is significantly lower than the value 0.935, sup- 
porting the presence of a sizeable frustration. Besides, 
with a smaller ratio, LGVO is expected to be closer to 
the fully frustrated point J2/J1 = 1/2 than LSVO, in 
qualitative agreement with the fact that no phase transi- 
tion was found in that system down to 1.9 K (see below). 
What we cannot say however on the basis of this analysis 
is which of the couplings J\ and J2 is larger. 

To be more quantitative, we need to know the ratio 
Tmax/© for the J\ — J2 model as a function of J2/ Ji- 
lt turned out to be impossible to get accurate estimates 
in the strongly frustrated region J2/J1 — 1/2, but ex- 
act diagonalizations with 3 sizes available (4, 8 and 16 
sites) give a reliable estimate for J2/J1 < .4, while Quan- 
tum Monte Carlo simulations, which suffer from the mi- 
nus sign problem, provide useful information down to 
J2/J1 — 2 on the other side, where only two sizes can be 
used for exact diagonalizations due to the type of order 
(8 and 16 sites). The experimental ratio T max /0 = .72 
for LSVO then implies that J2/ J\ is approximately equal 
to .1 or to 3.5, while for LGVO, T max /0 — .67 implies 
that J2/J1 is cither close to .25 or to 2.5. 

A significant difference between the two compounds is 
discernable if one reports the derivative dx/dT vs. T (see 
the inset to Fig. 2). One observes that around T c ~ 2.83 
K a peak is present for LSVO, while no anomaly in dxl dT 
is detected for LGVO, down to 1.9 K. The peak occurs 
at the same temperature where a peak in 7 Li NMR 1/Ti 
is observed (see Fig. 3), signaling a phase transition to 
a magnetically ordered state. Remarkably, T c was found 
independent on the magnetic field intensity, within ±0.15 
K (i.e. ±5% ), up to H = 7 Tesla. 

In LSVO, for H = 1.8 Tesla, one observes that 7 Li 1/Ti. 
is constant between 3.5 and 293 K. In the high temper- 
ature limit (T> 0), by resorting to the usual Gaussian 
form for the spin correlation function one has p0[ 



(VTi)c 



7 2 S(S +1)V2tt 



uj e 



k,i,j 



k \2 
ij) 



(1) 



with Aij (i = x,y,z, j = x,y) the components of the 
hyperfine tensor due to the k th V 4+ , 7 the gyromagnetic 
ratio and ui e — Jksy/2zS(S + l)/3/7i. z = 8 is the 
number of nearest neighbour spins of a V 4+ coupled via 
an effective superexchange coupling J related to J\ and 
J2. 7 Li hyperfine coupling constants in LSVO have been 
estimated by reporting the temperature dependence of 
the paramagnetic shift, both for the single crystal and for 
the powders, as a function of the susceptibility. It turned 
out that 7 Li nuclei are coupled to V 4+ ions both via a 
dipolar and a transferred hyperfine coupling At = 1700 
Gauss |ll|, which is attributed to the two V 4+ nearest 
neighbours. The dipolar term is close to the one esti- 
mated on the basis of lattice sums. From Eq. 1, by using 
the high temperature value of 1/Ti, one finds J = 6.4 
K, a value consistent with the Curie- Weiss temperature 
derived from the analysis of the susceptibility. 

For T < 0, V 4+ spins are strongly correlated and the 
behaviour of 1/Ti depends on which regime LSVO is: 
quantum critical, renormalized classical or quantum dis- 
ordered fl^ ]. However, one has to notice that the tem- 
perature dependence of 1/Ti can be strongly influenced 
by the g-dependent hyperfine form factor [ p^[ . Hence, 
in order to make significant statements on the correlated 
spin dynamics of this frustrated 2DQHAF we have first 
calculated the hyperfine form factor, which was found 
only weakly q— dependent in view of the sizeable trans- 
ferred coupling. Therefore, the temperature dependence 
of 1/Ti is fully determined by the one of the correlation 
length and the fact that 1 jT\ is constant down to ~ 3 K 
suggests that LSVO, for H = 1.8 Tesla, is in the quantum 
critical regime |l2| . 

Below T c one observes an activated decrease of 1/Ti 
which is typical of a magnetically ordered system with 
a gap in the spin wave spectrum [14] < An estimate of 
the gap amplitude can be done using a simple Arrhenius 
fit, which yields A = 18 K (see the inset to Fig. 3). 
Although the estimate of A might not be very accurate 
in view of a certain proximity to the phase transition it 
should be noticed that the estimated value is consider- 
ably larger than what one would expect for a magnetic 
system with a coupling constant of a few degrees kelvin, 
as deduced from susceptibility and 1/Ti measurements 
for T 3> 0. This fact could suggest a modification of 
the superexchange coupling at low temperatures, possi- 
bly involving a lattice distortion. The occurrence of a 
lattice distortion is, in fact, corroborated by the modifi- 
cations in 29 Si NMR powder spectra below ~ 3.4 K (see 
Fig. 4a). One observes, just above T c the appearence 
of a shifted narrow peak in 29 Si NMR powder spectrum. 
On decreasing temperature the low-frequency peak pro- 
gressively disappears while the intensity of the high fre- 
quency one increases. Two aspects should be remarked: 
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1) 29 Si NMR line does not broaden below T c indicating 
that the local field at 29 Si nuclei is zero; 2) the modifica- 
tion in the shift has to be associated with a modification 
in the chemical shift or hyperfine coupling, suggesting 
the occurrence of a structural distortion. The absence of 
a magnetic field at 29 Si site can be accounted for either 
by an AF order with the spins parallel to the c axis or 
by a collinear order with V 4+ spins along the sides of 
the square lattice. In order to exclude the presence of an 
AF order we studied the angular dependence of the mag- 
netic field at 7 Li nuclei below T c [|ll| . By considering the 
same hyperfine coupling tensor determined above T c 
we found that for the AF order the magnetic field inten- 
sity should always increase on turning the magnetic field 
from parallel to the c-axis to parallel to the ab plane, 
at variance with the experimental findings. Moreover, 
the fact that the spins are parallel to the ab plane below 
T c is supported by a recent EPR analysis of the g ten- 
sor []l6f , showing that there is a larger in-plane magnetic 
anisotropy. Therefore, wc conclude that the magnetic 
order is collinear with the spins in the ab plane. This 
is the first evidence of a collinear order in a frustrated 
2DQHAF, whose existence has been theoretically put for- 
ward long ago by Chandra and coworkers |^]/. It should 
be observed that only one of the two possble collinear or- 
ders, with a = ±1, is compatible with zero magnetic field 
at 29 Si, the one with the spins parallel to the staggered 
modulation, i.e. for spins along the x-axis the one with 
magnetic vector Q = (tt/cl, 0). The fact that in LSVO 
always one type of collinear order develops indicates that 
the four-fold symmetry of the square lattice was broken, 
possibly by the lattice distortion occurring just above T c . 

Further information on the collinear phase in LSVO 
can be derived from the temperature dependence of 7 Li 
NMR spectra below T c (see Fig. 4c). The temperature 
dependence of the order parameter, proportional to V + 
average magnetic moment, is obtained from the split- 
ting of the satellites of 7 Li NMR spectrum. The cen- 
tral peak and the two satellites do not correspond to the 
1/2 — * —1/2 and ±3/2 — > ±1/2 transitions, respectively, 
but correspond to Li sites where the local field is cither 
zero or non-zero (parallel or antiparallel to the external 
one). We have ruled out the possibility of a quadrupolar 
splitting by checking that both the length of the RF pulse 
yielding a it/ 2 rotation and the recovery curve of nuclear 
magnetization were the same for all lines. Moreover the 
observed splitting is nearly an order of magnitude larger 
than the quadrupolar one calculated on the basis of a 
point charge approximation. One notices (see Fig. 4b) 
a rather sharp, but continuous, decrease of the order pa- 
rameter close to T c . An accurate determination of the 
critical exponent [3 would require a temperature stabil- 
ity better than 5 x 10~ 3 K which could not be achieved 
with our cryogenic apparatus. Still an upper limit for f3 
can be estimated from our measurements which are car- 
ried out in steps of 1CT 2 K for T -> T c , We find that 



(3 < 0.25, a value compatible with a 2D Ising phase tran- 
sition, where (3 = 1/8. It should be observed that the 
relative amplitude of the central and satellite lines varies 
with decreasing temperature. This could be due to a 
modification of the interplanar correlation, to which 7 Li 
spectrum is sensitive. To elucidate this aspect further 
investigation of the collinear order is demanded. 

In conclusion, we have presented for the first time sus- 
ceptibility measurements showing that Li2VOSi04 and 
Li2VOGe04 can be considered as prototypes of frus- 
trated 2DQHAF and NMR spectra demonstrating that 
in the former a collinear phase is established at low- 
temperature, as predicted for J2/J1 > 1/2 ||. Finally 
29 Si NMR spectra suggest the occurrence of a structural 
distortion, just above the magnetic transition, which lifts 
the degeneracy between the two collinear ground states. 
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FIG. 1. a) Schematic phase diagram of a frustrated 
2DQHAF on a square lattice as a function of the ratio J2/J1 
of the superexchange couplings, b) Structure of LSVO and 
LGVO projected along [001]. (Si,Ge)04 tetrahedra are in 
grey, VO5 pyramids are in black while the grey circles indi- 
cate Li + position. For details see Ref. 8. 

FIG. 2. Temperature dependence of the magnetic suscepti- 
bility x = M/H, for H = 0.3 Tesla, in LSVO (closed squares) 
and LGVO (open circles) powders. In the inset the corre- 
sponding temperature dependence of the derivative dx/dT is 
reported. 



FIG. 3. Temperature dependence of 7 Li NMR 1/Ti for 
H — 1.8 Tesla along the c axis in a LSVO single crystal. 
In the inset we report the corresponding Arrhenius plot of 
the experimental data for T < 2.66 K. The error bars when 
not reported are within the circles. 



FIG. 4. a) 29 Si NMR powder spectrum in LSVO for 
H = 1.8 Tesla. The dotted lines mark the position of the 
peak at high and at low temperatures, b) Temperature de- 
pendence of the splitting between the low-frequency and the 
high frequency peak for 7 Li NMR spectra in c). The solid 
line indicates the behaviour expected for an order parameter 
with a critical exponent (3 = 0.25. c) 7 Li NMR spectra for 
H — 1.8 Tesla along the c axis of a LSVO single crystal, in 
the proximity of T c = 2.83 ± 0.005 K. 
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Fig.l (R. Melzi et al.) 
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